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Abstract 

We investigate homogeneous and isotropic cosmological models in scalar-tensor 
theories of gravity where two scalar fields are nonminimally coupled to the geom- 
etry. Exact solutions are found, by Noether symmetries, depending on the form 
of couplings and self-interaction potentials. An interesting feature is that we deal 
with the Brans-Dicke field and the inflaton on the same ground since both are 
nonminimally coupled and not distinguished a priori as in earlier models. This 
fact allows to improve dynamics to get successful extended inflationary scenarios. 
Double inflationary solutions are also discussed. 
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1 Introduction 

Extended gravity theories have recently assumed a prominent role in theoretical physics 
investigations since any unification scheme, as supergravity or superstrings, in the weak 
energy limit, or viable early universe cosmological models, as extended inflation, seem to 
base on them. 

Besides, any effective theory, where quantum fields are taken into account in a curved 
space-time, results in nonminimal couplings between geometry and matter (scalar) fields 
IJ. Furthermore, notwithstanding the fact that Einstein's general relativity is experi- 
mentally tested with high degree of accuracy, from solar system tests to binary pulsar 
observational data, it has become a peremptory necessity to consider alternative theories 
of gravity The issues is now: What kind of theory? The plethora of them is overwhelm- 
ing from higher-order gravity theories, to Kaluza-Klein multidimensional theories, to 
induced-gravity theories, to gauge theories with torsion. Several of them seem to be 
consistent with some quantum gravity effect in the weak energy limit but, till now, no 
one can be universally considered the "full" quantum gravity theory. 

Despite of this shortcoming, most of them have remarkable physical implications from 
cosmology to particle physics. A particularly relevant role is played in inflationary cos- 
mology where, from the early Starobinsky model to the more recent hyperextended 
models, non-standard theories have been widely used. In fact, the goal of every infla- 
tionary model is to generate a brief period in which the scalar factor of the universe, a(t), 
increases superluminally, i.e. a(t) > t. If a(t) grows by e 60 or more during this period, the 
horizon, flatness and monopole problems can be resolved. In addition, inflation generates 
energy density fluctuations which may be seeds for large scale structure formation. This 
features can be obtained in several alternative theories of gravity. 

However, designing a detailed microphysical model that accomplishes all of these goals 
has proven to be extremely difficult and several times the extended theories have to be 
adjusted to reach some partial goal. 

As a general scheme, one needs a grand entrance into inflation which is a mechanism 
which drives the universe into a false vacuum phase. The large, positive vacuum density 
acts as an effective cosmological constant which triggers a period of a (quasi) de Sitter 
expansion. Then one needs a graceful exit: a mechanism capable of terminating the in- 
flationary expansion, reheat the universe to a high temperature, and restore a Friedman 
expansion. This second issue can result highly problematic due to the fine-tuning of pa- 
rameters which one needs to connect in the isotropy of the cosmic microwave background 
(e.g. the so called "big bubble problem" 

Cosmological models deduced from nonminimally coupled theories of gravity (e.g. 
Brans-Dicke or induced gravity) have provided schemes capable of escaping such difficul- 
ties. Extended inflationary models could accomplish several goals which earlier models 
failed and, in particular, they cure some shortcomings of "old" inflation. Extended infla- 
tionary models also start when the universe is trapped in a false vacuum state by a large 
energy barrier during a first-order phase transition ||. The failure of old inflation was 



1 



that the universe could never escape the false vacuum state since the rate of tunneling 
through the barrier remains small compared to the inflationary expansion rate ||. Ex- 
tended models, in the various versions, avoid the same failure by introducing mechanisms 
so that the tunneling rate exceeds the expansion rate and, hence, the transition to the 
true vacuum can be completed. For a comprehensive review, see ||. 

The key ingredient on which extended inflation lies is the relation between the tun- 
neling or bubble nucleation rate, A, and the expansion rate (the Hubble parameter) H. 
The ratio is the dimensionless "bubble nucleation rate" e = A/if 4 . The false vacuum can 
be percolated by true vacuum bubbles only if e exceeds a critical value, e cr i t — 0.2. In 
old inflation, e is time-independent since A and H (which depends on the false vacuum 
energy) do not vary during inflation [J/J. 

Two situations are possible: i) e < e crit , in which case the true vacuum bubbles 
never percolate and the universe inflates forever, or ii) e > e cr i t , in which case the true 
vacuum percolates, but so quickly that there is insufficient inflation to solve cosmological 
problems. A way to bypass this shortcoming is to avoid bubble nucleation altogether. 
New inflation and chaotic inflation utilize this approach. For example, in new inflation, 
the energy barrier disappears altogether as the universe supercools and the universe 
evolves slowly but continuously from the false to true vacuum phase || 10 1 . However, 
the model has to be fine-tuned. 

Extended inflation models employ an alternative approach asking for the time varia- 
tion of e. Initially, e is much less than e CT u to achieve sufficient inflation, but then it grows 
during inflation to a value e > e cr n so that the phase transition can be completed. Being 
A fixed by the form of the self interacting potential, the only quantity on which one can 
act is H . Cosmological solution with a time-varying H can be easily obtained by using 
modified theories of gravity like scalar-tensor theories. A simple extended inflationary 
model Q can be contructed by using the Brans-Dicke theory ||11|| . The gravitational 
action is 
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where <fi is the Brans-Dicke scalar field, uj is a dimensionless parameter and C m is the 
matter Lagrangian density including all the non-gravitational fields (from now on, we 
shall use natural units h = c = kg = 8ttG = 1). The cosmological Friedman-Robertson- 
Walker (FRW) equations are 

} + (1.3) 

where H = a/a, p m and p m are, respectively, the energy and pressure densities of matter, 
and finally, k = 0, ±1. For p m = —p m and k = we have 

^) = (l + f) 2 , (1-4) 
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a(t) = [l + ^) , (1.5) 



a 

where % i s an integration constant connected to the energy of false vacuum, a = (3 + 
2cj)(5 + 6cj)/12. 

Immediately we see that for yt < a, <p ~ 0, we have a de Sitter solution. If, for 
example, to > 90, one obtains the 60 e-foldings necessary to solve cosmological problems 
of standard model. When yt > a {t) evolves as a power-law expansion and H ~ t . 
This feature allows the successful graceful exit since e > e cr u. Of course, the main 
ingredients are the variation of Newton constant and the coupling of the geometry to the 
scalar field. In other words, the model succeeds because the effective Newton constant 
G e ff = _1 is decreasing, then H is decreasing and e is increasing. 

The main flaw of this model is related to the expected value of the parameter u. 



In order to restore Einstein's general relativity, we should have uo — > oo [12], then the 



value of u) in constrained by the classical tests of general relativity: light deflection and 



time-delay experiments require ui > 500 E| while the bounds on the anisotropy of the 



microwave background radiation give u < 30 fll3 |. In conclusion, to must be a function 
of time in order to obtain viable models. A pure Brans-Dicke theory is not able to yield 
realistic models and we have to introduce, at least, a function u = uj(4>) in order to 
overcome the above difficulties. Several proposal have been done to improve the early 
extended inflationary model and, in some of them also A is assumed to vary 0. The 
strict condition which implies A =constant is a feature connected to Brans-Dicke models 
as it is shown in ||. 

In the so called hyperextended inflation [15| 



one can assume 



L0(<f>) = Uo + u; m( p m , (1.6) 



to improve the Brans-Dicke model 22 . If m = 5, the microwave background bounds are 



satisfied flEfl . Alternatively, assuming the coupling 

*= F ^- wW = mm? ' (L7) 

one can get successful implementations without fine-tuning the initial conditions. In 
particular, if F(ip) is a sixth order polynomial, the big bubble problem is avoided since 
the model is independent of the bubble-size distribution |L6] . 



Other approaches give interesting results. For example, it is possible to include a 
first- or second-order potential V(4>) for the Brans-Dicke field as in induced gravity 



theories |[L7| . In this case, the potential places constraints on the percolation time-scale 
in the graceful exit and, furthermore, it can give rise to multiple episodes of inflation 
which may reveal extremely useful for large scale structure formation (e.g. super-cluster, 
cluster and galaxies). 

Another way to escape the u;-parameter constraints is to consider a curvature-coupled 



inflation |18||. Also in this case, extended inflation results enhanced since, in some sense, 
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the roles of Brans-Dicke field and inflaton are mixed. This feature allows to satisfy the 
solar system constraints on lu, to avoid the big-bubble problem, to construct models with 
double inflationary episodes. 

A more sophisticated way to bypass the graceful exit problem can be obtained by 
coupling first-order phase transitions to curvature-squared inflation J0J . The mechanism 



(getaway inflation) is based on a nonminimally coupled higher-order gravity theory where 
terms like ip 2 R 2 appear in the usual gravity-inflaton action. Their role is to produce an 
inflationary phase of the background which has a classical end. At the same time, a stage 
of bubble production via semiclassical tunneling occurs allowing useful spectra for large 
scale structure formation. 

A final remark concerns the role which the Brans-Dicke scalar could have for dark 
matter in extended inflation. Its oscillations in the various models could account for the 
discrepancy between the dynamical estimate of the density of matter in the universe, 
Q ~ 0.2, and the prediction of inflation, Q = 1 [ED] which seems to be confirmed by the 



BOOMERANG experiment [21 



All these arguments and several more make extremely interesting to search for cosmo- 



logical solutions useful for extended inflation. A first investigation in this sense is in p2 | 
where general scalar-tensor theories of gravitation were studied in order to "model" use- 
ful extended inflationary behaviours. Intermediate inflationary universes with expansion 
scale factor of the form 

a(t) = a expt p , 0<p<l, (1.8) 

were found. These models allow to succeed in realizing phase transition and graceful 
exit. 

More recently, Modak and Kamilya |23[ derived exact cosmological solutions by the 



so called Noether Symmetry Approach [pj] in scalar-tensor gravity theories discussing the 
role of the coupling function u(ip) connected to the Noether symmetry. They improved 
the approach in ||25|| , where symmetries and solutions were found for theories with a scalar 
field nonminimally coupled to gravity, by introducing a second scalar field (the inflaton) 
as in extended inflationary models. Exponentially expanding solutions, in asymptotic 
region, were found and this feature does not allow to solve the graceful exit problem also 
if general relativity was asymptotically recovered. 

In this paper we want to discuss, by Noether Symmetry Approach, a further gen- 
eralization taking into account two nonminimally coupled scalar fields and their self 
interaction potentials. In this way, the roles of the Brans-Dicke field and the inflaton are 
mixed and both fields are taken on the same ground. This fact could be coherent with 
the stochastic approach for the fundamental laws of nature since the role of the fields is 
not attributed a priori p6 | . 



The paper is organized as follows. In Sect. 2, we discuss the double scalar-tensor 
action, derive the equations of motion, the point-like FRW Lagrangian and the cosmo- 
logical equations. Sect. 3 is devoted to the Noether Symmetry Approach which has to 
be improved for the double field case since the configuration space results enlarged. The 
summary of found symmetries is given considering also the subcases where one and not 
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two nonminimal couplings are present. The cosmological solutions are given in Sect. 4 
while the graceful exit problem is discussed in Sect. 5. Conclusions are drawn in Sect. 6. 



2 Double Scalar- Tensor Action and Equations of Motion 

The most general action in four dimensions, where gravity is nonminimally coupled to 
two scalar fields noninteracting between them, is 



F(<p)R + G(^)R + - V(<p) + ±W - W(iJ>) 



(2.1) 



where we have not specified the four functions F(ip), V((p), G(if)), and W(ip)- This action 
generalizes those used till now to construct extended inflationary models^. The Brans- 
Dicke action (|1 . 1| ) can be immediately recovered by using the transformations (|1.7|) . In 
our units, the standard Newton coupling is recovered in the limit F(ip) + G(ip) — > —1/2. 
The field equations can be derived by varying with respect to g^ u 

[F(<p) + G{n (R pv - l -g, u R) = T$> + T#> . (2.2) 

In the right hand side of (|2.2|) there is the effective stress-energy tensor containing the 
nonminimal coupling contributions, the kinetic terms and the potentials of the scalar 
fields ip and if), that is 

T$> = ~<P^<P-,u + \g^aV' a ~ \g^V{ip) - g^UF{ip) + F(<p) ;ta , (2.3) 
and analogously, 

T$ = ~M ;u + \g^;a^' a - \g^W{^) - g^Gty) + G{^ v . (2.4) 

□ is the d'Alambertian operator. The variation with respect to ip and ip gives the Klein- 
Gordon equations 

(dF\ dV , , 

and 

, n fdG\ dW , , 



1 We point out that the more general action is 

A = J d^xj—g [F(cp, t(,)R - V{(p, 1>) + A(<p, V-)(V^) 2 + B((p, ^)(V^) 2 
but for the purpose of the paper, we will confine ourselves to the action (|2.l[) 
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Their sum is equivalent to the contracted Bianchi identities |25 |. Let us now take into 
account a FRW metric of the form 



ds = dt —a 



dr 2 



kr' x 



(2.7) 



and substitute it into the action ( |2.1| ). Integrating by parts and eliminating the boundary 
terms, we get the point-like Lagrangian 



dF - 3 '- 2 
£ = 6——a 2 dip+6Fad 2 —6kFa 

dip 



arip 



^^-a 3 V(ip)+6—a 2 dip+6Gad 2 -6kGa+ — 
2 dip 2 



3„/,2 



The Euler-Lagrange equations, corresponding to the cosmological Einstein equations are 

2 k 



+ 



[F + G] 
,d 2 F dF 



a a 
2-+ - 

a \a 



+ 



a 

dF 



+ 2 



.dF +rt pdG 
^ dip dip 



+ 



-a 3 W{ip). 

(2.8) 
is are 

(2.9) 



dip 2 

2 

+ 6 



| - 2 d 2 F | ~ 

^ dip dip 2 dip 



.(^ + ^ 2 )-(V + W) 



o, 



' .dF + ^dG 
^ dip dip 



- ) + ^[F + G] + -(ip 2 + ip 2 ) + V + W = , (2.10) 



ip + 3[-)ip + 6 



ip + 3 [ - ) ip + G 



d\ 2 k 

-) + ~ 

a) a z 



a fa 
- + (- 

a \a 



+ 



'dF\ dV _ 

K dip J dip 

'dG\ | dW _ Q 

^dip J dip 



(2.11) 



(2.12) 



where Eq. ( |2.10| ) is the energy constraint corresponding to the (0, 0)-Einstein equation. 
Let us now go to solve the system (|2.9|) - (|2.12|) by using the Noether Symmetry Approach. 
The solutions strictly depend on the form of the functions F, G, V and W. By the Noether 
symmetries it is possible to select these functions so that the system ( |2.10| )-( p.l2j ) can 
be reduced and then integrated. 



o Selecting Couplings and Potentials by the Noether Symmetries 

Given an undefined extended gravity theory, the existence of a Noether symmetry can 
select the form of the coupling and the scalar field potential [^J or the form of the 
higher-order Lagrangian density, e.g. f(R, OR) J27j. 

At the same time, as we will show in the next section, the symmetry allows to reduce 
the dynamical system by a cyclic variable making it easier to solve. Taking into account 



the Lagrangian Q2.8D , its configuration space is three-dimensional, Q = {a,ip,tp}. In the 
language of quantum cosmology, it can be identified with a minisuperspace |[25|| . The 
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tangent space on which the Lagrangian Q2.8| ) is defined is TQ = {a, a, ip, ip, ip, ip} so that 
the lift vector X, the infinitesimal generator of symmetry, is 

X ^_|_^^_|_ d ^ da d ^ df3 d ^ d 
da dip dip dt da dt dip dt dip ' 

where a, (3, 7 are functions of a, ip, ip. A Noether symmetry exists if the condition 

L x C = 0, (3.2) 



is realized. Lx is the Lie derivative wit respect to X. Properly speaking, Eq. ( |3.2j ) 
corresponds to the contraction of the vector X with the Lagrangian ( |2.8|) . The constant 
of motion connected to the Noether symmetry is nothing else but 

S = i X #c (3.3) 

where 

dC , dC , dC , , /n AS 

$c = 1 rrda+—dp+- T diP 3.4 
da dip dip 

is the Cartan one-form given by a Lagrangian £ and ix is the contraction with respect to 
X. The relation between Eqs. ( [3.2|) and Q3.3|) can be easily seen if the vector is generally 
expressed as 

x = Q <A + !|!s (3 . 5) 

dq l dt do 1 



Using the Euler-Lagrange equations, it can be shown that 25 



If the Noether symmetry exists, Eq. ( |3.6|) gives (|3.2| ). In the Hamiltonian formalism, by 
a Legendre transformation, we get 

H = air a + ipTTp + tp-ir^ - C , (3.7) 

where TT q = dC/dq, q = {a, ip, ip} are the conjugate momenta. The phase-space vector 
for the symmetry is now 

• 9 . d ■ d .. d .. d yd , n oN 

1 = a d~a + Vty + a % + a d~d + ^ + ^ ' (3 ' 8) 

and a Noether symmetry exists if 

L r H = . (3.9) 
The conserved quantity (|3.3|) and the Hamiltonian (p.7|) gives the Poisson brackets 



{S ,H} = 0. (3.10) 
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Our issue is now to determine the functions F, G, V, W by this Noether symmetry tech- 
nique. We shall adopt the Lagrangian formalism. The condition ( |3.2|) gives the system 
of partial differential equations 



F\a + 2a 




3a + 12 



3a + 12 



d 2 F 

a/3 + 



dip 2 
d 2 G 



07 



d^ 



+ 



' da dp 
2a + a— + a— 
^ da dp 

' 9a; 07 

2a + a— — h a— 
, aa cfyy 
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.da 
dp 




dF\ da dp 

dp J dp dp 

dG\ da ^ $7 

dip J dip dip 



da „ a 2 9/5 07 

+ 2 -7T- ^ + — TT + a TT" 

099 6 aa <9<^ 

^ (9a a 2 07 9/5 
ch/> 6 9a 90 




dip) dp 




(3.13) 



0, (3.14) 



0, (3.15) 



da 
dip 



'dF\ dp 07 
v aV / dip dp 







(3.16) 



6A; 



(dG\ 
a + a& + 7 — a 



+a z 



3aV + /3a + 3alU + 7a 



(3.17) 

obtained by equating to zero the second degree coefficients in a, p, ip. The number of 
these equations is 1 +n(n+ l)/2, where n is the dimension of the configuration space Q. 
The system ( |3.11| )- (|3.17| ) is the straightforward generalization of the system (5.30)-(5.34) 
in |25| for the case of two scalar fields. 

The integration of ( |3.11| )-( |3.17| ) gives as a result the functions a(a,p,ip), P(a,p,ip), 
•y(a,p,ip), F(p), V(p), G(ip), W(ip). Solutions are not unique and the various cases 
which we have found are summarized in Table I. 

In Table II, being G(ip) = 0, the cases where only one nonminimal coupling is present 
are summarized. The value of the spatial curvature constant k is also given. The quanti- 
ties F , Go, Fq, Ai j2 are constants. Using these results, the dynamical system ( |2.9|) -( pT2"D 
can be reduced since, as we shall show below, a change of variables can be found where a 
cyclic coordinate is present. This feature allows to integrate more simply the dynamics. 



4 The Cosmological Solutions 

The existence of a Noether symmetry gives, in any case, a cyclic variable so that the 
transformation 

C(a, a, p, p, ip, ip) * — > C{w, w, u, u, z) (4.1) 
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is always possible. If more than one symmetry exists, more than one cyclic variable 
can be present f25f . In a geometric language, it is always possible to choose a new set 
<?i = Qi(Qk) i,k = 1, 2, 3, adapted to the foliation given by X 

i x dQ 3 = l, i x dQ 3 =0, j = 1,2, (4.2) 

where i x , as before, is the contraction given by X and dQj = (dQj/dqi)dqi. Explicitly, 
in our case, Eqs. ( |4.2j ) become 

dw „dw dw 

a ^ + %^w- =0 - (4 ' 3) 

dz n dz dz , 
a Ya + % + ^ = 1 ' (4 ' 5) 

where z is the cyclic variable. However, the transformation (O)-(O) are specified as 
soon as the functions a, (3, 7 are given. As an example, let us take into account Case 1 
in Table I. Cases 2 and 3 of Table I and 1,2,3 of Table II can be deduced from it. The 
system ( |4.3j )-( |4"75| ) is solved by the choice of the new variables 

w = a 3 ip 2 , u = c?ip 2 , z = In a. (4.6) 

For the scalar-field functions, as we said, we choose Case 1 in Table I. Lagrangian (|2.8| ) 
becomes 

1 1 

C = 6F z(w-2zw) + — (w-3zw) 2 -V w + QG z(u-2zu) + — (u-3zu) 2 -W u (4.7) 

8w Sw 

while the equations of motion are 



U ZiXL J 

(11) W ^ \ 

--7T^)+8Vo = 0, (4.9) 
w 2w z J 

(8F - l)w - (32F - 3)zw + (8G - l)u - (32G - 3)zu = S (4.10) 

Q{8Fo-l)zw-3{32Fo-3)z 2 w+—+8V w+6{8Go-l)zu-3{32G -3)z 2 u+ —+8W u = , 

w u 

(4.11) 

where, clearly, z is the cyclic variable and So is the constant of motion connected to 
z. With respect to the system ( [2.9|) -( fTT2j ), system (|4.8| ) - (|4.11| ) is reduced and it is 
highly symmetric due to the functions F((p), G(ip), V(<p), W(ip) selected by the Noether 
symmetry. 
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Since the role of the two fields is completely symmetric, we can suppose that, de- 
pending on the value of the parameters, one can select two regimes of physical interest. 
For example, in the first case, dynamics is <p (-u)-dominated, in the second case, it is 
ip (w)-dominated. The relation among the initial data is given by the energy condition 
( |4.11| ). We get in the (p (w)-dominated regime 



z(t) = Z\ arctan 
for Cic 2 > 0, and 



u(t) = ci exp[A]t] + c 2 exp[— Ait] 



■ exp [Ait] + z 2 In |ci exp [Ait] + c 2 exp [-Ait] | + z 



(4.12) 
(4.13) 



,Ci 



z(t) = — ^arctanh yJ\ — | exp[— Ait] 1 + z 2 In \c\ exp[Ait] + c 2 exp[— Ait]| + z , (4.14) 



for cic 2 < 0, while in the ip (w)-dominated regime 

w(t) = c 3 exp[A 2 t] + c 4 exp[-A 2 t] , 

z(t) = z 3 arctan ^^^exp[A 2 t]^ + z 4 In |c 3 exp[A 2 t] + c 4 exp[— A 2 t] | + z 
for C3C4 > 0, and 



(4.15) 
(4.16) 



, C 3, 



z(t) = — ^ 3 arctanh yJ\ — |exp[— A 2 t]J + z 4 In |c3 exp[A 2 t] + C4 exp[— A 2 t]| + zq , (4.17) 



for C3C4 < 0. The constants Zq, Ci, Zj, i, j = 1, . . . , 4 are integration constants while 



Ai 



^(32^0 - 3) 
4F (12F - 1) 



A, 



W (32G - 3) 
4G (12G - 1) 



(4.18) 



assume the role of cosmological constants depending on the couplings and the potentials. 
Inverting the relations fl4.6|) , we get 



Z\ arctan 



a(t) = ao (ci exp[Ait] + c 2 exp[— Ait]) exp 

tp(t) = (f (ciexp[Ait] + c 2 exp[-Ait]) -2 ^exp 
for Cic 2 > 0, and 

a(t) = ao (ci exp[Ait] + c 2 exp[— A 2 t]) exp ziarctanh 



Z\ arctan ( J — exp [Ait] 



'— exp[Ait] 

C2 



(4.19) 



, (4.20) 



\ — I 

c 2 



exp [Ait] 



(4.21) 
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ip(t) = ifQ (ci exp[Ait] + c 2 exp[-A 2 t]) 2 ( exp 



(4.22) 



for C\C2 < in the ^-dominated regime. 

The situation is analogous in the ^-dominate regime, but the constants Ai, Zi, ci )2 
and ipo have to be substituted with A 2 , z 2 , c 3i4 and ip . 

It is easy to see hat we have two inflationary eras. Their durations are ruled by 
the parameter Ai,A 2 which strictly depends on the strength of the couplings. Another 
interesting particular solution is 

a(t) = aoe*', ^) = ^oexp{^^ I t} , rj,{t) = ^ exp {g^ 2 ^*} , (4-23) 



Vo = - 0/ J l ^ F (f -±)( Fo --±-) (4.24) 



where 

L ' ~ ~8(8F - Lf u V^ u 12 J V" 10.6 

^ = -8^8ctT) G »( G «-^)( G »-T^)- < 4 ' 25 > 

Also here the inflationary behaviour is clear. 

By using similar arguments, we can analyze Case 2 in Tab. I Here the potential terms 
cancel each other in the Lagrangian fl£7|) . We get the power-law solution 

a(t)=a t\ tp{t) = ip Q r n , m = tfj r n , (4.26) 

which is particularly useful for extended inflation being n an arbitrary constant depending 
on the initial conditions. A similar situation holds in Case 4, which is a minimally coupled 
case with two fields and two cosmological constants. For k = and Ai = — |A| 2 , one 
finds 

a < ( > = °° ( "- ^ = ^<rhnj tl ' 3 " +Vo ' = (4 ' 27) 

for n 7^ 1/3, and 

a(t) = a t 1/3 , cp(t) = -j\nt + (f , ^(t) = lnt + , (4-28) 

when n = 1/3. In this case, two Noether symmetries are present and they assign the 
value of gravitational constant being 

Fo + Go = _5M+M). (4 . 29) 

Let us now analyse, in detail, Case 5. Without loosing of generality, we can assume 
Fq = F = 0, 70 = 1 and studying the couplings F((p) = (p 2 /12 and G(ip) = G . 
Lagrangian ( |2.8|) becomes 

£ = + 6G ^j aa 2 + a 2 a^ + a 3 - lij + (4.30) 

11 



where A = Ai + A 2 . Clearly if) is the cyclic variable. Eqs. (|4.3|) - (|PD are satisfied by 

w = a , u = cup — if) , z = cup . 

With the further change 



X = z -u, 



Lagrangian ( |4.30|) reads 



C 



6G w 2 + 



z 

y 



w + 



X 



-A r 



where two cyclic variables appear. The dynamical system is 

•2 



2 



x 2 



QG (2ww + w 2 ) = — + 3w 2 y - A 

Xw 3 = Si , 
zw = S 2 , 
z 2 + 12G w 2 + x 2 w 2 + 2Aw 2 = , 
whose general solution is given by the elliptic integral 

w 2 dw 



vMiw 6 + A 2 w 2 + A 3 



±t 



where 



A, 



A, 



A, 



Y 2 

^2 



-4, 



6G 6G 12G 

In the particular case where A 2 = 0, we get the explicit solution 



a(t) = a 
<p(t) = 



sinh(±3J^M) 



^ sinh(±3v^Tt) 



cosh(±3v/]4T) i 



3V]47 sinh^iv^M) 



In tanh 



^±3v^M N 



(4.31) 
(4.32) 



(4.33) 



(4.34) 

(4.35) 
(4.36) 
(4.37) 



(4.38) 
(4.39) 

(4.40) 
(4.41) 

(4.42) 



which asymptotically gives a de Sitter behaviour. A last interesting case is 8 in Tab. I, 
which can be assigned by the functions 



F = F , V{<p)=A : 



G(if>) = G + ^- + 
lb 



k = -1 

W(ifj) = 



1P 4 
Wo 



A 



(4.43) 



being G{if>) and W(if>) free for the Noether symmetry. The model is relevant for hyper- 
extended inflation (see e.g. ii~5fl). Power law solutions like in [|l|, [U| are easily found. 
The cases in Tab. II are essentially subcases of those discussed above. 
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5 Inflation and graceful exit 



As we discussed in Introduction, the goal to get a sufficient inflationary period and then 
to exit from it, without imposing any particular fine tuning, can be achieved by assuming 
the variation of the bubble nucleation rate e. However, we are taking into account a first 
order phase transition after which we recover a Friedman stage ||. In principle, being 
e = X/H A , we can expect the variation of both A and H. The form of A strictly depends 
on the form of the theory and, as it is discussed in ||, it is time independent toward 
the late times, if we are dealing with a Brans-Dicke theory. In our cases, by using ( |1 . 7|) , 
we get that most of the couplings selected by the existence of Noether symmetry can be 
recast in a Brans-Dicke-like form. 

In spite of the variation of the effective gravitational coupling, it is reasonable to 
assume A to be approximatively constant |J, so that the mechanism of the graceful exit 
can be essentially connected to the variation of H. 

However, this argument does not work for more general classes of theories, as hyper- 
extended inflation |15| , where w(4>) is not a constant. Among the cases in Tabs. I and 



II, we have also couplings of the form 

see e.g. the cases 5 in Tab. I and 4 in Tab. II. This situation deserves more attention 
since we can distinguish a regime where we match a sort of hyperextended inflation 
(<p — ► 0) and a regime where the extended inflationary scheme is recovered (<p — > oo). 
In any case, the microwave background bounds have to be satisfied, as discussed in [jT6| . 
Furthermore, taking into account double-field models, the contributions to A come from tp 
and ip. Being both fields nonminimally coupled, and from the forms of couplings selected 
by the Noether symmetry, we are dealing with a double Brans-Dicke-like theory where 
the extended inflationary mechanism is improved. Looking at the solutions of previous 
section, we can have double inflationary stages ruled by the parameters of couplings and 
self-interaction potentials (see e.g. ( |4.19| )-( fL22| ) or ( |4.26| )). This situation is extremely 



interesting since "very" large scale structure and large scale structure can be selected 
by these inflationary phases. In fact, we can have "two" first-order phase transitions 
and then "two" bubble nucleations where the size of bubbles is given by the coupling 
parameters. In other words, we can expect two graceful exits given by the superposition 
of two extended inflationary phases. To be more precise, at a given time t > to after 
nucleation, the "comoving" bubble radius is 

r(t,t ) = f dta{t)~ l , (5.2) 

while the "physical" size of the bubble is 

K(t,t ) = a(t)r(t,t ) . (5.3) 
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When t — > oo, the form of a(t) selects the size of the bubble. In the cases ( |4.19| ),( fL22| ) 



and (f4.23|) , this size is finite since the asymptotic behaviour is a(t) ~ exp H t. For power 



law behaviours, the growth of the bubble size is linear. 

Besides, we have a variation of the Hubble parameter in most of the cases we dealt 
with: in ( |419| ), ( |4.21| ), and ( }4.40| ), it converges to a constant for t — > oo, in ( |4.23| ), it is 
exactly a constant, in the other cases, it is H ~ t -1 . 

Graceful exit is achieved if, being A a constant, e is less than e crit during inflation and, 
after bubble nucleation, e > e cr #. In our cases, H is the key parameter which governs the 
behaviour of e. For power-law solutions, as for standard extended inflationary models, 
the graceful exit is easily recovered (see Eqs( |4.26| ) and ( |4.27| )). In the asymptotically 



exponential cases, the parameter e goes to a constant for t — > oo and the graceful exit 
problem has no solution. In fact, the function H, calculated from ( |4.19| ) and ( |4.21 



is 



a sort of step function with two different constant values at t — > ±00. The related H~ A 
has a singularity in the origin which does not allow a graceful exit from inflation. The 
situation for the solution ( |4.40|) is similar. 



6 Conclusions 

In this paper, we derived exact cosmological solutions in double scalar-tensor gravity 
theories by the general approach of searching for Noether symmetries. This work gen- 



eralizes those in |23], |24], Pq] . The couplings and the potentials of both scalar fields are 
connected with the existence of the symmetries, and the solutions of dynamics furnish 
power law or de Sitter evolutions. As a consequence, in all the above cases it is easy to 
calculate the bubble nucleation rate e = A/if 4 to test if one succeeds in graceful exit. 
Depending on the value of intervening parameters, this can be accomplished in several 
cases. 

Furthermore, being in principle, both fields nonminimally coupled and self-interacting 
with a potential, their role is mixed and it is not possible to distinguish, a priori, a Brans- 
Dicke field and an inflaton field as in other extended inflationary models. This distinction 
seems, in our opinion, rather artificial in view of a stochastic approach to the fundamental 
interactions where the effective role of the various fields is distinguishable only in the low 



energy limit (see [g(J and reference therein) and there is no reason why a field should 
interact with the gravitational field and the other one not. 

Another remark deserves double inflation which is ruled by the parameters of the 
theory and, then, by the Noether symmetry. As it is well known this feature is of 
extreme interest in perturbation theory since it can furnish the seeds for the formation 
of structures at large and at very large scales. As we have seen, it is very common in our 
approach and it could contribute to the enhancement of a successful extended inflation. 

Finally, we want to stress the fact that the standard Newtonian coupling can be 
recovered in several of the above models, being 
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so that as soon as F((p) — > F and/or G(ip) — > Go, general relativity is restored (in our 
units Fq + Go — > —1/2) and both fields can contribute to its recovering. This means 
that in an accurate setting of the models, one could succeed both in graceful exit and in 
recovering of standard gravity. 



Table I - Symmetries in double nonminimally coupled models. 



N. 


a 


P 


7 


F(tp) 


G(^) 1%) W(^) 


fc 


1 


a 


-3y?/2 


-3^/2 


F <p 2 


GoV 2 Vfa 2 W„^ 2 





2 


a 


-3v?/2 


-3ip/2 


F <p* 


G ^ 2 


A 


-A 





3 


a 


-3v?/2 


-3^/2 


FoV 2 


G ^ 2 





W 2 





4 








7o 


F 


Go 


Ai 


A 2 


VA; 


5 





1/a 


7o 




Go 


Ai 


A 2 





6 





1/a 


7o 


±<p* + Fq(p + F 


G 





A 





7 





1/a 


7o 




Go 


Ax 


A 2 





8 





A) 





F 


G(V) 


0,A 




VA; 


9 





1/a 





±<f? + FqV + F 




0,A 








Table II - Symmetries in single nonminimally coupled models. 



N. 


a 


P 


7 








1%) 




fc 


1 


a 


-3/2y? 


-3^/2 




F cp 2 











2 


a 


-3v?/2 


-3^/2 




F y? 2 












3 


a 


-3y?/2 


-3^/2 




FoV 2 




A 


-A 





4 





1/a 


7o 


t^ 2 









A 





5 





1/a 


7o 


tV 




- i*o 


Ax 


A 2 





6 








7o 








Ai 


A 2 


VA; 
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